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Abstract 

We use a modified SU(2) chiral sigma model to study nuclear matter component 
and simple bag model for quark matter constituting a neutron star. We also study the 
phase transition of nuclear matter to quark matter with the mixed phase characterized 
by two conserved charges in the interior of highly dense neutron stars. Stable solutions 
of Tolman-Oppenheimer-Volkoff equations representing hybrid stars are obtained with 
a maximum mass of 1. 67 Mq and radius around 8.9 km. 
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1 Introduction 



The existence of quark matter in neutron stars have been considered by many authors ' 
over the past two decades. In a simple but fairly accurate representation, a neutron 
star is constituted mostly of neutrons and a small number of protons whose charge is 
balanced by the leptons. In works prior to 1992, local charge neutrality was imposed 
on all possible phase transitions above saturation density, leaving the systems with only 
one conserved quantity, namely the baryon number. Such a system is described as 
one component system. However, in 1992, Glendenning^ pointed out that neutron star 
matter has two independent conserved components, namely the baryon number and the 
electric charge and hence charge neutrality must be applied as a global and not a local 
condition. For a one component system, the transition takes place at a constant pressure 
where the two phases having different densities, for all proportions of the two phases are 
in equilibrium. The equilibrium pressure is found through Maxwell construction^. For 
a two-component system, the common pressure and all properties of the two phases in 
equilibrium vary with the proportion of the two phases and this region of two co-existing 
phases in equilibrium is usually referred to as the mixed phase^. 

In this paper, we have extended our earlier one component analysis^ of neutron star 
structure and associated phase transition using a two component formalism. We consider 
the baryon chemical potential (/xb) and electric charge chemical potcntial(yU£;) as the 
two conserved charges. Using Gibb's criteria for a two-component system a first order 
phase transition between beta-stable nuclear matter and quark matter is indicated where 
pressure and density vary continuously in the mixed phase region. Taking the existence 
of such a phase transition between nuclear matter and quark matter as a guide, we solve 
the Tolman-Oppenheimer-Volkoff(TOV) equations with appropriate nuclear matter^ and 
quark matter^ equations of state and find the hybrid stars to consist of a quark-matter 
core with the nuclear matter forming a thin crust. 

This paper is organized as follows. In sec. 2, we present the equation of state(EOS) 
for nuclear matter. Then we proceed for equation of state of quark matter presented in 

1 



sec. 3. In sec. 4, phase transition is discussed. In sec. 5, we discuss the structure of hybrid 
stars. Finally we discuss and summarize our results in sec. 6. 



2 Hadronic Phase 

We have used the modified SU(2) chiral sigma model(MCH)^ for hadronic phase since 
chiral model has been very successful in describing high density nuclear matter. The 
importance of chiral symmetry* in the study of nuclear matter was first emphasized by 
Lee and Wick^ and has become over the years one of most useful tools to study high 
density nuclear matter at the microscopic level. The nonlinear terms in the chiral sigma 
model give rise to the three-body forces which become significant in the high density 
regime^''. Further, the energy per nucleon at saturation needed the introduction of 
isoscalar field^^ in addition to the scalar field of pions^^. We also include the interaction 
due to isospin triplet p- vector meson to describe the neutron rich- matter 

The modified SU(2) chiral sigma model'' considered by us includes two extra higher 
order scalar field interaction terms which ensures an appropriate incompressibility of 
symmetric nuclear matter at saturation density. Further, the equation of state(EOS) 
derived from this model is compatible with that inferred from recent heavy-ion collision 
data^l 

The EOS for hadronic phase is calculated by using the Lagrangian density^, 
AC 



8m4 

+l9'yu;,u;^ - \g,..G^^ + \mlp,.p^ - \gMP.-ri')i> (1) 

In the above Lagrangian, F^^ = d^uj^ - d^uj^, (7^^ = - d^Pfj, and x = (tt^ -Fcr^)^/^, -0 
is the nucleon isospin doublet, tt is the pseudoscalar-isovector pion field, a is the scalar 
field and a;^, is a dynamically generated isoscalar vector field, which couples to the 
conserved baryonic current = '07/i'0- is the isotriplet vector meson field with mass 



rUp. B and C are constant coefficients associated with the higher order self- interactions 
of the scalar field . 

As obtained in our earlier work^, the total energy density(e) and pressure(P), for the 
/9-stable nuclear matter is given by 



^ m^l - y-f ^ _ C _ 



The energy per nucleon is ^ = ^ and the chemical potential \s ji = {P + e)/nB- The 
concentrations of neutrons, protons and electrons can be determined using conditions of 
beta equilibrium and electrical charge neutrality^ i.e /in — IJ'p + IJ'e — Ue , (here /li 
is the chemical potential of the particle species i). 

The values of five parameters Co-, C^^, Cp, B and C occurring in the above equa- 
tions are obtained by fitting with the saturation values of binding energy /nucleon (-16.3 
MeV), the saturation density (0.153 fm~^), the symmetric energy(32 MeV), the effec- 
tive(Landau) mass (0.85M)^^, and nuclear incompressibility (~300 MeV), in accordance 
with recent heavy-ion collision data^^ are C^, = 1.999 fm^, = 6.816 fm^, Cp = 4.661 
fm^, B = -99.985 and C = -132.246. The equation of state for the present model is found 
to be softer^ with respect to the original one^^ 



3 Quark Phase 

We employ the bag model in our study of the quark phase. We consider here the quark 
matter EOS which includes u, d, and s quark degrees of freedom in addition to electrons. 
We have taken the electron, up and down quark masses to be zero^ and the strange quark 
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mass is taken to be 180 MeV^'^^. As shown in Ref. 5, The energy density and total 
pressure of quarks including the bag constant(B) is given by 

P^Pe + Y.Pf-B (3) 

/ 

e = -^ + Ew- (4) 

/ 

Where is the electron pressure and Pf is the pressure for quark flavor with / = u, 
d or s. The chemical equilibrium condition i^d — IJ's — IJ'u + A*e written^ in terms of 
baryon(//B) and electric charge(/XE) chemical potential renders P and e as functions of 
two independent chemical potentials fiB and fiE- Defining the number density of quark 
flavor f as Hf = the baryon number density and charge density of the quark system 
is given by 

riB = Jl^^/ 
/ 

P = 3(2n„ -rid- ris) 
Charge neutrality of this system requires 

3^*^ — ^ 

4 Phase Transition 

For the phase transition from /^-stable nuclear to quark matter, we consider the above 
described hadronic phase and quark phase with a typical value of ag — 0.5, B^/^ = 150 
MeV^'^^. The Maxwell construction within one component formalism was done^, by as- 
suming that the two phases are separately charge neutral keeping pressure constant with 
discontinuity in energy and number density over the transition region. The transition 
was found to extend over the density range of 0.27 - 0.37 fm~^. 

In the present work we employ theoretically more accurate two-component analysis 
which makes the pressure and density vary continuously over the mixed phase region. 
We follow Glendenning procedure^ of parameterising nuclear matter by two chemical 
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potentials {ij>b,I^e) i-e baryon and electric charge chemical potentials. Gibb's condition 
for phase equilibrium at zero temperature for both phases gives 

PnpifJ'B, fJ'E) = PgpifJ'B, IJ'e) (5) 

As described in Ref. 17, the nuclear and quark matter share a common electron distribu- 
tion in the mixture is more appropriate. Hence the overall charge neutrality requirements 
is that the nuclear, quark and lepton components together are charge neutral, i.e 

XPQpil^B, I^e) + (1 - x)Php(i^b, I^e) = Peil^B, I^e) (6) 

Where x is the filling fraction of quark matter in the mixed phase. The above two 
equations(5,6) involve three variables /ib, IJ'E and So for a given value of /ib, the two 
equations can be solved simultaneously to give the values of and x- Then the energy 
density emp and baryon density hmp of the mixed phase can be calculated using the 
equations^^ 

cmp = X^QP +{^-x)^HP 
riMP = XTT-QP + (1 - x)nHP (7) 

This procedure is carried out provided that < x < 1. For x = the system is in pure 
hadronic phase and for x = 1 the system is in pure quark phase. 

For as = 0.5, B = (150 MeV)'', it is found that the number density and energy density 
of the mixed phase extends over the range 0.20 - 0.58 fm~^ and 191 - 638 MeV/fm~^ 
respectively. The total pressure thus obtained is now no longer constant over the density 
range of the mixture; it increases monotonically in the mixed phase region as shown in 
Fig.l by solid line. In the current picture of bulk phases, the local charge neutrality of 
both phases is relaxed to a global charge neutrality condition which indicates that in 
the mixed phase region the nuclear matter phase will be positively charged whereas the 
quark matter phase will be negatively charged^. As the phase transition implies that the 
neutron star consists of quark matter in the interior and nuclear matter in the periphery 
hence it can be called as a hybrid star. 
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Figure 1: Pressure (P) as function of baryon number density in the mixed phase region. 
The solid curve corresponds to two-component system and the dotted for one component 
system. 

5 Structure of Hybrid Stars. 

As a hybrid star is assumed to be a spherically symmetric distribution of mass in hydro- 
static equilibrium, the equilibrium configurations are obtained by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equations^^ for the pressure P and the enclosed mass m, 



dP{r) 
dr 



-G 



[e(r) + P{r)\[m{r) + 47rr^P(r)] 
r^[l — 2Gm{r)/r] 



(8) 



dm{r) 
dr 



47rr^e(r) 



(9) 



Where G is the gravitational constant, m(r) is the mass contained in a volume of radius 
r. In order to get stellar radius R and the gravitational mass M = m{R) we integrate 



the TOV equations for a given central energy density e{r — 0) — Cc with the boundary 
condition m(r = 0) = and P{r = i?) ~ 0. 

We have plotted in Fig. 2 the mass of hybrid star as a function of central energy 
density to examine the stability of such a star. Taking into account the stability of such 
stars under density fluctuations require dM/dCc > 0^°. As may be seen from the figure, 
dM/dec becomes negative around 1800 MeV/fm^ after which it may collapse into black 
holes^^'^°. This yields the maximum mass of hybrid star as M ~ 1.67Mq. Fig. 3 shows 
the mass as a function of radius obtained for different central densities for such a star 
which indicates the maximum radius to be around 8.9 km. 

We observe that for a given central energy density of 700 MeV/fm^ such a star has 
pure quark matter core of radius 3.29 km, mixed phase covers 5.08 km and pure nuclear 
matter crust of 0.21 km, whereas for ec = 1800 MeV/fm^, the corresponding datas are 
6.69 km, 1.63 km and 0.1 km. Hence it is clear that if we take smaller central energy 
density then pure quark phase is found to cover less region and the mixed phase covers 
more region out of the total radius of the hybrid star and the reverse happens in case 
of larger central energy density. But pure nuclear matter covers a thin crust of the 
hybrid star in the whole range of relevant central energy densities. Thus the existence 
of a broad mixed phase and a much thinner layer of pure nuclear matter phase (~ |th 
of the one-component system) are the two most important features distinguishing the 
characteristics of a two-component system vis-a-vis a one-component one. We have 
exhibited the results of both analyses in all our figures for quick visual comparision. 

Since the data with regard to symmetric energy(32±6 MeV) and nuclear incompress- 
ibility (300±50 MeV) which have been used to fix the coupling constants of the theory 
are associated with some errors, we have also tried to check the sensitivity of our results 
for Mq and R with respect to these errors. To do so we have computed the coupling 
constants (C^, Cu,, C(t,B and C) for four cases of symmetric energy value 32 MeV with 
incompressibihty 250 and 350 MeV and the incompressibility of 300 MeV with symmet- 
ric values 26 and 38 MeV. The maximum values for M and R are found to be M = 
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Figure 2: The mass(M/MQ) of the hybrid star as a function of central energy 
density (cc). The solid curve corresponds to two-component system and the dotted for one 
component system. 

1.67^0 oo^M© and R = 8.9lo;o3 kms. Such small effect of errors in symmetric energy and 
incompressibility on maximum values of M and R can be attributed to the fact that only 
about 2% of the radius of the hybrid star contains pure nuclear matter. 

We also calculate the surface gravitational redshift Zg of photons which is given by^^ 

Zs = , ^ - 1 (10) 

^(1 - 2GM/R) 

In Fig. 4, we have plotted Zg as a function of M/Mq . In this context it may be mentioned 
here that our result for the surface redshifts lying in the range of 0.2 to 0.5 agrees quite 
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Figure 3: Mass as a function of radius for the hybrid star. The solid curve corresponds 
to two-component system and the dotted for one component system. 

well with the values determined from gamma ray bursters^^. We then compute the 
relativistic Keplerian angular velocity given by^^ 



0.72, 







(11) 



as for neutron stars. Fig. 5 shows our result for variation of relativistic Keplerian angular 
velocity as a function of M/Mq for such a star. We find that the VLk has an inverse 
relationship after an initial increase with the mass of the star beyond 1.67M0. This 
indicates that we can not have mass of hybrid star more than about 1.67Mq. We also 
observe that the maximum value of VL^ being near about 1.18 xlO"^ sec~^ implies the 
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Figure 4: T/ie surface gravitational redshift(Zs) as a function of star mass (M/Mq). The 
solid curve corresponds to two-component system and the dotted for one component sys- 
tem. 

maximum period to be 0.53 millisecond that agrees with the results obtained by Burgio 
et.aP. 

6 Conclusions 

In this work we have extended our earlier one-component analysis^ of phase transition 
and structure associated with the hybrid stars using a two-component formalism. As 
expected, the results of both analysis share similar qualitative features with certain 
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Figure 5: The Keplerian angular velocity{D,k) as a function of star mass (M/Mq). The 
solid curve corresponds to two-component system and the dotted for one component sys- 
tem. 

quantitative differences. We find that a first order phase transition exists between the 
hadronic phase and quark phase at density of about four times the nuclear matter density. 
The phase transition from nuclear matter to quark matter indicates that the core of a 
neutron star consists of quark matter. To obtain stable hybrid star solution, we have 
solved TOV equations using appropriate equations of state and have taken as = 0.5 
and B=(150 MeV)^ for quark matter EOS. We observe that a stable hybrid star with 
quark core and a nuclear crust exists upto €c — 1800 MeV/fm^ beyond which instabihty 
is indicated. For the values of the parameters used in this model we found that the 
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maximum mass and radius of hybrid star to be about 1.67lo!ooi-^0 8-9iao3 km 
respectively. It is also observed that if we take smaller central energy density then pure 
quark phase is expected to cover less region and the mixed phase covers more region out 
of the total radius of the hybrid star and the reverse happens in case of larger central 
energy density. But pure nuclear matter covers a thin crust of the hybrid star in all 
situations. 

Using this modified SU(2) chiral model the maximum mass and radius of neutron 
star was observed in Ref.7 as M = 2.1Mq, R = 12.1 km. The corresponding results for 
hybrid star characterized by a mixed phase with one conserved charge obtained in our 
earher work^ was M = 1.69lo:ooiM0, R = 9.3lo:o5 km. Thus we find that a hybrid star 
characterized by a mixed phase is more compact than a neutron star. Further we notice 
that sharpening of the analysis of the structure of hybrid star with inclusion of two 
conserved charges instead of one, stabihses the star in a shghtly more compact form(M 
= 1.67io;oo[M0, R = 8.9lo;o3 km). Wc also find that two-component analysis leads to 
much larger mixed phase region inside the hybrid star with corresponding shrinkage of 
the pure nuclear matter crust. These we consider as important distinguishing features 
of the results obtained with two and one component formalism. Overall it is observed 
that there is a decrease in maximum mass due to the existence of mixed phase in the 
core of a hybrid star. The greater compactness of the star also leads to smaller time 
period lying in the sub-millisecond range as obtained by us. It is also observed that the 
surface gravitational redshift and relativistic Keplerian angular velocity of the hybrid 
star can not increase beyond M/Mq — 1.67; showing a decrease with increase in M/Mq 
beyond this value. 
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